Future supernova probes of quintessence 
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We investigate the potential of a future supernovae data set, as might be obtained by the proposed 
SNAP satellite, to discriminate between two possible explanations for the observed dimming of the 
high redshift type IA supernovae, namely either (i) a cosmological evolution for which the expansion 
of the universe has been accelerating for a substantial range of redshifts z ~ 1; or (ii) an unexpected 
supernova luminosity evolution over such a redshift range. By evaluating Bayes factors we show 
that within the context of spatially flat model universes with a dark energy the future SNAP 
data set should be able to discriminate these two possibilities. Our calculations assume particular 
cosmological models with a quintessence field in the form of a dynamical pseudo Nambu-Goldstone 
boson (PNGB), and a simple empirical model of the evolution of peak luminosities of the supernovae 
sources which has been recently discussed in the literature. We also show that the fiducial SNAP 
data set, simulated with the assumption of no source evolution, is able to discriminate the PNGB 
model from a number of other spatially flat quintessence models which have been widely studied 
in the literature, namely those with inverse power-law, simple exponential and double-exponential 
potentials. 



I. INTRODUCTION 

Type la supernovae (SNe la) can be used as standard 
candles to infer the luminosity distance, dr,, as a function 
of redshift 0-0] , and such data provide a key element in 
the case for cosmic acceleration. The measurements pro- 
vide a simple way to estimate cosmological parameters. 
This is the aim of at least two collaborations: the Super- 
nova Cosmology Project |l]|| and the High-redshift Su- 
pernova Search Team J2|. However, the analysis of some 
~ 50 supernovae at redshifts around z~0.5 has not yet 
provided strong constraints on the nature of the "dark 
energy" . 



"Quintessence" j4Tp| has been proposed as a candi- 
date for the "dark energy" to provide a dynamical so- 
lution to the cosmological constant problem R. The 
supernova data, although sufficient to constrain the pa- 
rameters of quintessence models IqHlQl , is not yet abun- 
dant enough to allow a discrimination between compet- 
ing quintessence models. Moreover, the unexpected dim- 
ming of the type IA supernovae at redshifts z ~ 0.5 is a 
phenomenon which could be readily attributed either to 
an accelerating Universe or to an unexpected luminosity 
evolution pl[ . The present data alone is not sufficient 
to allow a discrimination between these two possibilities 
with complete confidence [|l(4[l2||l5J 1 . 

Very recently the identification of a SN la event at red- 
shift z~1.7 p3 has provided tantalizing evidence that 
the expansion of the universe could have been deceler- 
ating at that epoch. Naturally, no firm conclusions can 
be drawn from a single event. However, if this result 
holds up it would rule out the simplest models of source 



evolution. Furthermore, supernovae events at such red- 
shifts are the type of observations which should enable a 
discrimination to be made between some of the different 
quintessence models. 

Much more data is needed to enable an accurate es- 
timation of the nature of the dark energy. This might 
be accomplished by a dedicated space telescope, the Su- 
perNova Acceleration Probe (SNAP) [|15|Jl6[ , which aims 
to collect a large number of supernovae with z < 2. In 
this paper, we assess the ability of the SNAP mission to 
determine various properties of the "dark energy" . By 
analysing a simulated data set, as might be obtained by 
the proposed SNAP satellite, we can test the ability of 
such experiments to distinguish among currently attrac- 
tive quintessence models p|p]j9| Jirj| , ^7| -p3| . 

The feasibility of determining the properties of the 
dark energy component by using simulated data sets has 
already been considered by several authors p4,p4 p8| . 
One common approach is to assume that the quintessence 
field is described by a perfect fluid with equation of state 
P = wp, where w is approximately constant over epochs 
of interest, or else slowly varying. For example, various 
authors J25|-p8fl consider models with an equation of state 
linear in redshift, w(z) — w$ + w\z. However, many real- 
istic cosmological models could fall outside the confines 
of these approximations. In this paper we wish to take 
a different approach by considering models in which the 
quintessence field is directly given by a Lagrangian, and 
in which w(z) has the freedom to vary widely over mea- 
surable redshifts. We will fit simulated SNAP data sets to 
the exact dr,(z) of different quintessence models derived 
from numerical integration of the coupled Einstein-scalar 



field equations. 

For the fiducial cosmological model used in the sim- 
ulation of the supernova data set, we consider a form 
of quintessence, an ultra- light pscudo Nambu-Goldstone 
boson (PNGB) |17j which is still relaxing to its vacuum 
state. Our reasons for this choice of quintessence model 
are twofold. Firstly, from the viewpoint of quantum field 
theory PNGB models are the simplest way to have nat- 
urally ultra-low mass, spin-0 particles and hence per- 
haps the most natural candidate for a presently-existing 
minimally-coupled scalar field. The effective potential of 
a PNGB field 6 can be taken to be of the form fl8|] 



II. 



ACCELERATING UNIVERSE VS 
LUMINOSITY EVOLUTION 



y(</»)=M 4 [cos(</»//) + l], 



(1) 



where the constant term ensures that the vacuum en- 
ergy vanishes at the minimum of the potential. This 
potential is characterized by two mass scales, a purely 
spontaneous symmetry breaking scale / and an explicit 
symmetry breaking scale M . 

The second motivation for choosing a PNGB model 
rather than other forms of quintessence is that it pro- 
vides a natural framework for studying the question of 
the possibility of source evolution. 

Analyses of the SN la data performed with source 
evolution have been undertaken in the case of open 
Friedmann-Robertson- Walker (FRW) models [jl3| . How- 
ever, the measurement of the positions of the first 
acoustic peaks in the angular power spectrum of cos- 
mic microwave background radiation anisotropies by the 
BOOMERANG, MAXIMA and DASI experiments now 
gives unequivocal evidence that the Universe is close to 
being spatially flat |2S|-|3l|]. This makes the choice of a 
spatially open FRW model uncompclling. 

On the other hand, the PNGB cosmologies have the 
virtue that while they are spatially flat, there is no a 
priori preference for either accelerated or decelerated ex- 
pansion [ |10|Jl2| . Both possibilities are available at mod- 
est redshifts, depending on parameter values. Ultimately, 
the scalar field will undergo coherent oscillations about 
its minimum, and the resulting luminosity distance will 
become indistinguishable from that of an Einstein-de Sit- 
ter model at late times, although the fraction of energy 
density in clumped matter, f2 m , and the fraction of en- 
ergy density in quintessence, £1$, can take any values 
consistent with Vt m + Cl$ = 1. Whether the scalar field 
is currently rolling down the potential for the first time 
in the history of the universe - leading to a luminosity 
distance relation similar to that produced by a cosmolog- 
ical constant - or whether it has already undergone one 
or more oscillations by the present epoch, is a matter of 
a choice of initial conditions. Beyond some bounds set 
by primordial nucleosynthesis these initial conditions are 
not very much constrained by our present knowledge of 
the models, resulting in diverse possibilities for cosmolog- 
ical evolution. The current and final values of £l m and fi^ 
likewise depend on initial conditions for the scalar field, 
and on the parameters M and /. 



In an earlier paper |10|, we considered the observa- 
tional constraints arising from SNe la data and gravita- 
tional lensing data on cosmological models based on Ein- 
stein gravity minimally coupled to a scalar quintessence 
field with a PNGB potential (|l|) . In the case of the super- 
novae, we studied the constraints on the PNGB param- 
eters M and / from 60 supernovae data from the Super- 
nova Cosmology Project (hereafter P98) of Perlmutter 
et al. j|. We numerically evolved the coupled Einstein- 
scalar field equations of motion forward from the epoch 
of matter-radiation equality to obtain the luminosity dis- 
tance redshift relation, comparing the results with the 
data both with and without an allowance for the possi- 
bility of source evolution. 

In allowing for source evolution we are acknowledging 
that the peak luminosities of distant SNe la have been 
normalized according to empirical "Phillips relations" 
P2| - f34|l between observed peak luminosity and super- 
nova decay time, which have been found to be valid at low 
redshifts. Although we would hope that such relations re- 
main applicable at high redshifts, until the Phillips rela- 
tions can be modelled and accounted for physically some 
doubt will always remain about applying these relations 
at higher redshifts. 

A simple empirical model for possible source evolu- 
tion was employed in ref. p0| : following Drell, Loredo 
and Wassermann |l3| we added a term /?ln(l + z) to the 
distance modulus. This particular luminosity evolution 
function is simply chosen as an illustrative example, and 
is not singled out by any particular physical model. Nat- 
urally, one can criticize it on these grounds. However, 
the Phillips relations are also purely empirical, and the 
purpose of our study as with that of ref. |13| was simply 
to test the extent to which any form of source evolution 
was able to account for the observed data as compared 
to an accelerated expansion. 

For the purpo ses of analysing the data it was assumed 
in refs. |10|,|l2|,|l3| that the prior for the parameter j3 was a 
Gaussian with a mean /3o and standard deviation b. The 
best-fit value of /3o for the P98 data set was found to 
vary very slightly with initial conditions, taking a value 
A, = 0.414 for w t = <j>i/f = 1.5 and f3 = 0.435 for 
w.i = 0.2, where <pi is the value of the quintessence field 
at the time of onset of matter domination [z = 1100). In 
this paper we use the initial value Wi = 0.2. Our reason 
for doing this is that in general, and especially for larger 
values of u>i , there is a tension between the values of Q^q 
and Hoto favoured by current observations, with larger 
values of i?o^o often corresponding to values of fi^o rather 
greater than 2/3. It is still possible to find parameter 
values of M and / which give Q^o ~ 0.7 while yielding an 
age for the Universe consistent with the lower bound of 
the current observational range |35| , even for u>i = 1.5. 
However, one finds that the tension between the values 



of fi^o and H^t^ is more easily mitigated by choosing a 
lower value of Wi , for which the scalar field spends more 
of its early dynamical history higher up the potential hill. 
Confidence limits on the M, f parameter space obtained 
from the P98 data set, assuming no evolution, are shown 
in Fig. |. 

In refs. @,|§ we included all 60 data from the P98 
data set. In fact, in ref. |2[] Perlmutter et al. considered a 
subset of 54 SNe la which excluded 6 supernovae events: 
two that are the most significant outliers from the average 
light-curve width, two that have the largest differences 
between observed and expected magnitudes (or fluxes), 
and two that are likely to be reddened. It was discovered 
that excluding these six supernovae produced a more ro- 
bust fit of the cosmological parameters £l m o and J7ao • In 
refs. |l(j,[l2| we took a conservative approach and included 
all 60 SNe la, since we are investigating models with dif- 
ferent cosmological parameters from those fitted in ref. 
j| and we did not wish to prejudge matters. However, 
one may of course redo the analysis of refs. Jl(],[l^] for 
the "fit C" data set of Perlmutter et al. For the PNGB 
model we find that in the best-fit case this reduces the 
normalized y 2 parameter from 101.6, i.e., 1.75 per degree 
of freedom, for the full data set to 58.65, or 1.13 per de- 
gree of freedom, for the reduced "fit C" data set. For the 
models with source evolution the normalized \ 2 is simi- 
larly reduced from 101.1, or 1.74 per degree of freedom, 
to 58.21, or 1.12 per degree of freedom. Thus the "fit 
C" data provides a more robust fit in all cases we have 
studied, whatever the model assumptions. 

Since the results in ref. [ [To|]i1| were based on the full 
P98 data set, we have included results based on both the 
full data set and the reduced data set in Fig. El to com- 
pare the differences. Note that although the 2ct allowed 
area to the right disappears entirely for the reduced data 
set for w — 0.2, this is not the case for w — 1.5, and 
thus the discrepancy between the P98 data set and that 
of Riess et al. g] which was noted in Ref. [1(|] remains 
for the PNGB model whether or not one uses the full or 
reduced data set. 

Here we will take a somewhat different approach to 
ref. Jul in our treatment of the case with source evolu- 
tion. We will assume that the parameter (3 has a flat 
prior bounded over some range A/3 = 2.3, with limits 
corresponding to —0.6 < [3 < 1.7, and marginalize the 
likelihood function over this prior. The choice of the 
bounds on /3 in explained below in section IIB. 

Despite the differences from the approach used in Ref. 
|[10f the conclusions of the analysis are broadly similar, 
and are shown in Fig. |2j. In particular, two regions of 
parameter space appear to be singled out: region I (see 
Fig. |l|), which corresponds to cosmological models for 
which the scalar field is rolling down the potential (Q) for 
the first time (with <f> increasing) at the present epoch, 
and region II (see Fig. @) in which the scalar field is rolling 
down the potential for the second time (with <fi decreas- 
ing) at the present epoch. In the absence of any evolution 
the P98 data set favours Region I, whereas if evolution is 



allowed for then region II appears to be favoured slightly 
more than region I. 

Note that over the entire parameter space considered, 
we find the value of (3 with the best fit of all is (3 *=* 0.622 
using the entire P98 data set, which occurs at parameter 
values M = 4.55 x lO^/i^eV, / = 0.884 x 10 18 GeV, 
while for the "fit C" data set it is (3 sw 0.720, which 
occurs at M = 4.64 x lO^/i^eV, / = 0.901 x 10 18 GeV. 
Averaged over the (M, f) parameter space we find ((3) rs 
0.286 for the full P98 data set, and 0) sw 0.348 for the 
"fit C" data set. Thus in both cases the value of j3 is 
peaked in region II, but for the "fit C" data the overall 
values of (3 are somewhat larger. 

The results of refs. JlC],[l3] were inconclusive. Os- 
tensibly the source evolution models appeared to pro- 
vide a slightly better fit. However, the extent to which 
they were better was not quantified statistically. One 
relatively straightforward way to compare the relative 
strengths of rival theoretical models in fitting a com- 
mon set of data is by the calculation of Bayes factors 
[]13t , p6|| , even when the models involve different numbers 
of parameters. Given more parameters, one may always 
find a better fit to the data - however, the Bayes factor 
approach includes an effective "Ockham's razor factor" 
which adds a weighting against the inclusion of extra 
parameters. The Bayes factor method may be readily 
applied to the results of |l0|,|l2[ to give an explicit statis- 
tic to compare the relative strengths of a cosmological 
explanation of the SNe la data versus an explanation in 
terms of source evolution. 

In Bayesian inference, when comparing rival models 
Mi, each with parameters Qi, the likelihood for a model 
conditional on the data (D) in a model comparison cal- 
culation is equal to the average likelihood for its param- 
eters, 



p(D\Mi) = \ dQ % p(Qi\Mi)C{Qi 



(2) 



where p(Qi\Mi) is the prior probability for Qi, and £(Qi) 
is the sampling probability for D presuming Mi to be 
true. The ratio of model likelihoods, 



Bi 



p(D\Mj) 
p(D\M d ) 



(3) 



is called the Bayes factor. When the prior odds does not 
strongly favour one model over another, the Bayes factor 
can be interpreted just as one would interpret an odds in 
betting. Kass and Raftery [p6| provide a comprehensive 
review of Bayes factors, and the recommended interpre- 
tation is summarized in Table 1. 



A. P98 data set 

In order to ascertain what degree of improvement could 
be expected with data from SNAP, we will begin by cal- 
culating Bayes factors for the existing P98 data set. 
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FIG. 1. Confidence limits on M,f parameter values, with w% = 0.2: (a) for the 60 supernovae la in the P98 data set; (b) 
for the 54 supernovae la in the reduced "fit C" P98 data set. Parameter values excluded at the 95.4% level are darkly shaded, 
while those excluded at the 68.3% level are lightly shaded. Overplotted are the contours for fi^o (dashed) and Hoto (dotted). 
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FIG. 2. Confidence limits on M,f parameter values, with Wi = 0.2, marginalized over a flat prior for the luminosity 
evolution parameter /3: (a) for the 60 supernovae la in the P98 data set; (b) for the 54 supernovae la in the reduced "fit C" 
P98 data set. Parameter values excluded at the 95.4% level are darkly shaded, while those excluded at the 68.3% level are 
lightly shaded. Overplotted are the contours for n^o (dashed) and Hoto (dotted). 



The prior ranges for parameters play an important role 
in Bayesian model comparison. We will make a choice 
of the prior ranges of parameters in the different mod- 
els by the resulting values of f2 m o and i?o^o that they 
give rise to. In particular, we will choose a conserva- 
tive bound of 0.5 < f^o < 0.9. For the expansion 
age of the universe we choose 0.72 < Hota < 1.15: for 
Hq = 70 km s _1 Mpc _1 it corresponds to a conservative 
bound of t o — 13 ± 3 Gyr consistent with recent esti- 
mates |p5| . The parameter space bounded by these two 
constraints is indicated by the contours of f2 m o an d H t Q 
in Figs. [0 and 0. In the large / region in region I where 
the contours appear to be parallel to the /-axis, a cut-off 
at / ~ 10 19 GeV, the Planck scale, is chosen. 

We calculate the average parameter likelihood by inte- 
grating C(M, /) over the prior ranges for the parameter 
and dividing the integral by the relevant area. For the 
full P98 data set the Bayes factor for the model without 
luminosity evolution versus the model with luminosity 
evolution, which is the ratio of the average parameter 
likelihoods, is B ~ 2.6. For the reduced "fit C" data set 
we obtain a Bayes factor of 2.5, which is the same as far 
as its interpretation is concerned. Thus the data alone 
cannot discriminate between the two hypotheses. 



B. SNAP data sets 

We will now simulate a data set which would be ex- 
pected to be obtained by the SNAP satellite to inves- 
tigate its potential for discriminating between the two 
possibilities. We will simulate two data sets from fiducial 
cosmologies chosen with parameters in each of regions I 
and II. In each case we will introduce a random error to 
the exact distance moduli to simulate a future supernova 
data set that has been converted to a table of effective 
magnitudes mi and redshifts z; of objects with a sin- 
gle fiducial absolute magnitude Mo- We will consider 
both statistical and systematic uncertainties in the mag- 
nitudes. Typically the redshifts are known to sufficiently 
high precision that their uncertainties can be ignored. We 
assume the supernovae observed uniformly within four 
different redshift ranges with the following different sam- 
pling rates, which are the same as those chosen by Weller 
and Albrecht [M: In the first range from z = 0-0.2 we 
assume that there are 50 observations, in the second and 
largest redshift range from z = 0.2-1.2 there are 1800 
SNe and in the two high redshift bins, z — 1.2—1.4 and 

TABLE I. Interpretation of Bayes Factors 



Bij 


Strength of evidence for Hi over Hj 


1 to 3 


Not worth more than a bare mention 


3 to 20 


Positive 


20 to 150 


Strong 


> 150 


Very Strong 



z = 1.4-1.7, there are 50 SNe and 15 SNe observations re- 
spectively. The statistical error in magnitude is assumed 
to be c mag = 0.15, including both measurement error 
and any residual intrinsic dispersion after calibration. 

For data set A, we assume there is no luminosity evo- 
lution. We choose the fiducial parameters M — 2.53 x 
10" 3 /i 1 / 2 eV and / = 0.58 x 10 18 GeV; these parameters 
give (fl m0 , r^ ) w (0.23, 0.77) and H t « 0.9. For data 
set B, we assume there is a luminosity evolution and we 
add a term /? ln(l + z) to the distance moduli. We choose 
the fiducial parameters M = 4 x 10^ 3 h 1 ' 2 eV and / = 
0.676 x 10 18 GeV. These parameters give (f2 mO ,O0 O ) ~ 
(0.25, 0.75) and H t ~ 0.8. The best-fit value of the 
parameter (3 for this particular choice of M and / is 
$ = 0.659 for the full data set, and $ = 0.727 for the 
"fit C" reduced data set. We will take j3 ~ 0.659 as the 
fiducial parameter for data set B, as it is the somewhat 
more conservative value, representing less luminosity evo- 
lution. 

The fiducial parameters for the two data sets have 
been deliberately chosen so as to give similar values of 
f2 m o,f^o and Hat{) on one hand, so as to be consistent 
with values favoured by other observational tests. On 
the other hand, the fiducial parameters for data set A 
are centred in region I, which corresponds to a universe in 
which the scalar field is rolling down the potential for the 
first time at the present epoch, whereas the fiducial pa- 
rameters for data set B are centred in region II, which cor- 
responds to a universe in which the scalar field is rolling 
down the potential for the second time at the present 
epoch. As was discussed in [J10| region I is favoured if the 
SNe la luminosity distances are true cosmological dis- 
tances, whereas region II becomes significantly favoured 
if the data hides a simple luminosity evolution. 

We will use Bayesian analysis to estimate uncertain- 
ties. Rigorous calculation of the likelihood for the 
quintessential parameters Q is very complicated, requir- 
ing the introduction and estimation of many additional 
parameters, including parameters from the light-curve 
model and parameters for characteristics of the individ- 
ual supernovae. With several simplifying assumptions 
fll3|, the final result is relatively simple. One finds 



£(G) 



e -x7a 



where 



x 2 (Q) = E 



\fx s<i - ufa; Q)f 



fi,i 



In the above equation, 

ft(zi',Q) = 5logdL{zi]Q)+25 



(4) 



(5) 



(6) 



is the distance modulus predicted by each model with 
parameters Q, while fi Sy i = rrn — Mo is the simulated dis- 
tance modulus, and its uncertainty is a^i — 0.15. Note 

that we fix Hubble parameter to Hq — 70 km/s Mpc 
to simplify the computation. 



The model with an unexpected luminosity evolution 
corresponds to replacing eq. (|5|) with 

&q,0) = y: [mm " W1 + j i] ~ Kzf ' Q)? ■ (7) 



a 



fi.i 



We will marginalize over (3 with a flat prior, to obtain 
the marginal likelihood 



£(fi) = ^g/'/-'"- X '"' 2 



(8) 



The above integration can be performed by an analytic 
marginalization |10| , |l3| . We separate q from x 2 where 



«(G) 






E 



[h{z l -Q)-P(Q)\n{l + z l )f 



(9) 



is independent of /3. The integral is thus an integral over 
a Gaussian in (3 located at (3 with standard deviation a. /3 
is the best-fit value of (3 given Q, and a is its conditional 
uncertainty. They are given by 



-=y 

n-2 Z^ 



[ln(l + Zl )f 



KQ) = « 2 



E 



h(«i;fi)h(l + «i) 



where 



fe(«i; Q) = M«,i -Mte; 2) 



(10) 

(11) 

(12) 



As long as /3 is inside the prior range and ct <C A/3, the 
value of the integral is well approximated by a^/^jr, so 
that 



£(Q) = 



A/3 



-9/2 



(13) 



The question of what bounds to place on the prior 
range for the parameter /3 is an interesting one, since we 
are dealing with a purely empirical model of source evo- 
lution with no a priori restrictions on (3, other than that 
very large values of (3 would not be physically plausible. 
We will therefore seek the narrowest range of values of (3 
consistent with our numerical results. One condition on 
the prior range for (3 is that (3 also lies inside the prior 
range: setting too tight a bound on A/3 would count 
against parameter values of M and / which might oth- 
erwise be included. By explicit numerical integration we 
have found that for the initial conditions and the range 
of parameter values of M and / we have considered (3 
lies in the range —0.6 < (3 < 1.7. We will therefore take 
these bounds to be the prior range for /3. 



If we fit the data set obtained from fiducial model 
A assuming that there is no luminosity evolution, the 
parameters M and / are well constrained by the simu- 
lated data. The 95.4% confidence level contour bounds 
a very small region around the best-fit values of M sw 
2.55 x 10- 3 /i 1 / 2 eV and / « 0.592 x 10 18 GeV, with 
X 2 = 1910. These best-fit parameters are very close to 
the fiducial parameters that generate the data set. The 
average parameter likelihood is £ avo ~ 7.1 X 10 _4 £o, 
where Co — exp(— 1910/2). If luminosity evolution is 
considered to exist (see Fig. 0), fitting the data set from 
fiducial model A we find that the 95.4% confidence level 
contour includes a larger region in region I of the pa- 
rameter space around the fiducial parameters, and a re- 
gion in region II. The average parameter likelihood is 
£ ave ~ 1.9x 10~ 3 £o- Hence, the Bayes factor for the 
model with luminosity evolution versus the model with- 
out luminosity evolution is B ~ 2.7. We conclude that it 
is not possible to discriminate between the two hypothe- 
ses even though the underlying simulated data does not 
have any luminosity evolution. 

If we fit the data set obtained from the fiducial model 
B by assuming that there is no luminosity evolution, 
then the 95.4% level bounds a very small region around 
the best-fit values of M « 2.27 x 10" 3 /i 1 / 2 eV and / w 
0.764 x 10 18 GeV. The average parameter likelihood is 
£ avo ~ 9.4 x 10- 18 £ , where £ = exp(- 1924/2). If lu- 
minosity evolution is considered to exist (see Fig. Q), on 
the other hand, then the 95.4% level bounds a small re- 
gion around the best-fit values of M w 4.00 x 10~ 3 /i 1 / 2 eV 
and / w 0.695 x 10 18 GeV, with (3 w 0.626 and normal- 
ized x 2 = 1924. The average parameter likelihood is 
£ ave ~ 8.3 x 10~ 4 £o- This time, the Bayes factor for the 
model with luminosity evolution relative to the model 
without luminosity evolution is much greater than one, 
B ~ 9.1 x 10 13 . Thus if the true data were derived from 
fiducial model B, it would show up clearly as giving a 
(1l(z) which can be discriminated from the model with- 
out luminosity evolution. 

This result in no way contradicts the inconclusive re- 
sult obtained from data set A. If the underlying data 
were to truly follow a simple luminosity evolution, then 
the analysis above shows that this would stand out as a 
very strong positive signal in the Bayes factor test. If the 
underlying data has no luminosity evolution on the other 
hand, then it may be possible to obtain a good fit with 
an extra luminosity evolution parameter, simply because 
one has an extra parameter to fit. Thus a simple lumi- 
nosity evolution is easy to rule in, but difficult to rule 
out. Combined together the two results demonstrate the 
efficacy of the Bayes factor approach. With data set B we 
would obtain a decisive result, but by Ockham's razor we 
should reject the more complicated evolutionary hypoth- 
esis if we were to obtain an inconclusive result such as 
that pertaining to data set A. In fact, we have found that 
by narrowing the range of the prior A/3 that the Bayes 
factor for data set A can even be made slightly greater 
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FIG. 3. Confidence limits on M,f parameter values, with Wi = 0.2, marginalized over a flat prior for the luminosity 
evolution parameter f3, for the 1915 supernovae simulated assuming M = 2.53 x 10 _3 /i 1 ' 2 eV, / = 0.58 x 10 18 GeV, and (3 — 0. 
Parameter values excluded at the 95.4% level are darkly shaded, while those excluded at the 68.3% level are lightly shaded. 
Overplotted are the contours for Q.^ (dashed) and Hoto (dotted). 
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FIG. 4. Confidence limits on M,/ parameter values, with W; = 0.2, marginalized over a flat prior for the luminosity 
evolution parameter (3, for the 1915 supernovae simulated assuming M — Ax 10 -3 ft 1 ' 2 eV, / = 0.676 x 10 18 GeV, and /3 = 0.659. 
Parameter values excluded at the 95.4% level are darkly shaded, while those excluded at the 68.3% level are lightly shaded. 
Overplotted are the contours for fi^o (dashed) and Hoto (dotted). 



than the cut-off value of 3 listed in Table 1. Clearly any 
weakly positive results should therefore also be treated 
with some caution. 



III. PNGB VERSUS OTHER POTENTIALS 

The constraints on cosmological parameters which 
arise from type la supernovae have been broadly stud- 
ied for a wide range of quintessence models |^ |K||r2f . 
However, all these models seem to have ranges of pa- 
rameters that fit the existing data equally well, with no 
particular model standing out as being observationally 
preferred. The lack of data and the large statistical error 
bars are among the contributing reasons. In this section, 
we will investigate the potential of the SNAP supernovae 
data set to discriminate between competing quintessence 
models, as contrasted with the P98 data set, assuming 
no source evolution. 

We consider some simple scalar field potential func- 
tions which have been widely studied in the literature. 
They are the simple exponential potential P,PlflT 



V = V^e - *"* , 
the inverse power-law potential p|,pO[ 

V = v A( t>- a , 

and the double-exponential potential |9|]2 
V = V j4 cxp(-Ae v ^^) , 



(14) 



(15) 



(16) 



where n is related to Newton's constant by k 2 = 8ttG. 
We evaluate the luminosity distance - redshift relation 
by numerically evolving the Einstein-scalar field evolu- 
tion equations forward from the time of onset of matter 
domination, as what we did for the PNGB models. For 
the simple exponential and the double-exponential po- 
tentials, without lost of generality we choose the initial 
condition fa = 0. For the power-law potential, we begin 
the integration by assuming that the scalar field has al- 
ready reached the tracker solution. These leave us with 
two parameters in each case. In Figs, [g, g, and 0, we dis- 
play the confidence limits for the 54 supernovae la in the 
P98 "fit C" reduced data set, on the parameter spaces 
for these quintessence models. (The equivalent figures 
for the full data set have been omitted, as the la and 
2cr confidence regions are very similar to those shown in 
Figs. U"]?], but just slightly narrower.) 

In this section we will only compare models (|14|)-(|lq) 
on the basis of assuming that the true data follows the 
PNGB model, with the fiducial parameters of data set 
A. Naturally, we could assume any of the other poten- 
tials p4|)-(|l6|) as our fiducial model, and then compare 
the other models on the basis of a new fiducial data set. 
This would allow us, for example, to compare the inverse 
power-law potential ( |l5|) w ith the simple- and double- 
exponential potentials (|l4|), (ttq), a test which we have 
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FIG. 5. Confidence limits on the parameter space for the 
simple exponential potential model V — V"Ae _A ' t * ', for the 54 
supernovae la in the P98 "fit C" reduced data set. Parameter 
values excluded at the 95.4% level are darkly shaded, while 
those excluded at the 68.3% level are lightly shaded. Over- 
plotted are the contours for f2^o (solid) and Hoto (dotted). 




FIG. 6. Confidence limits on the parameter space for the 
inverse power-law potential model V = Va0 _c \ for the 54 
supernovae la in the P98 "fit C" reduced data set. Parameter 
values excluded at the 95.4% level are darkly shaded, while 
those excluded at the 68.3% level are lightly shaded. Over- 
plotted are the contours for Sl^o (solid) and Hoto (dotted). 
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FIG. 7. Confidence limits on the parameter space for the 
double-exponential potential model V = VAexp(— Ae ^), 
for the 54 supernovae la in the P98 "fit C" reduced data 
set. Parameter values excluded at the 95.4% level are darkly 
shaded, while those excluded at the 68.3% level are lightly 
shaded. Overplotted are the contours for fi^o (solid) and 
Hoto (dotted). The lower right-hand region is excluded from 
the plot due to computational difficulties. 



not performed here. However, we will restrict our at- 
tention to the fiducial data set A relevant to the PNGB 
potential (Ef) , since the amount of computer time involved 
in computing the fiducial data sets is large, and an anal- 
ysis based on a fiducial PNGB model will prove to be 
sufficient to illustrate the dramatically increased discrim- 
inatory power of a SNAP dataset as opposed to presently 
available SNe la data. Naturally, one could extend the 
discussion to other fiducial data sets, but we will leave 
that to future work. 



A. P98 data set 

Firstly we will calculate the Bayes factor of the P98 
data set to compare the PNGB model with all the other 
quintessence models. In order to calculate Bayes fac- 
tors which compare different theoretical models we must 
naturally make choices of the range of prior parameters 
integrated over, and we are dealing with different param- 
eters in the different models. Different choices of priors 
for these parameters would lead to some variation in the 
resulting Bayes factors. We will make a choice of the 
range of prior parameters by using the resulting values 
of O^o and i?o^o that they give rise to. However, the pa- 
rameter spaces are not completely bounded by the f^o 
and Hoto constraints. Without a physical cut-off for the 
parameters, such as a cut-off for / at Planck scale for the 



PNGB model, we will investigate in each case the depen- 
dence of the Bayes factor comparison on the prior ranges 
of the parameters. We will show that in some cases, the 
choice of the prior parameter ranges do not affect the 
conclusion from the Bayes factor analysis. 

For the simple exponential potential model, as shown 
in Fig. ||, the Q^o and Hoto contours diverge at large Va, 
putting a bound on the parameter space. In the small 
A region, where the potential is effectively a cosmologi- 
cal constant model whose value depends on Va only, the 
fJ^o and H to contours remain near parallel. For the full 
P98 data set the average parameter likelihood over the 
region bounded by the fi^o and Hoto constraints, with 
a cut-off at A = 0.1, is £ avo ~ 0.67£_r, where Cr is the 
average parameter likelihood for the PNGB model. In 
the small A region the likelihood only changes slightly 
with decreasing A and tends to a constant, correspond- 
ing to an average parameter likelihood which we estimate 
to be £ aV o — 1-5£r. Since the small A region extends to 
an infinite range in terms of the parameterization shown 
in Fig. pL if we enlarge the prior range of A to include 
smaller values, then it has the effect of bringing the aver- 
age parameter likelihood closer to the average parameter 
likelihood in the small A region, which is effectively an 
upper bound. Thus the Bayes factor for the simple expo- 
nential potential model versus the PNGB model lies in 
the range 0.67 < B < 1.5 for the full P98 data set. Sim- 
ilarly 0.65 < B < 1.5 for the reduced "fit C" data set. 
In each case the larger value is the one appropriate to 
including the entire small parameter range for A. Such 
a Bayes factor is too small to confidently discriminate 
between the two models. 

For the double-exponential potential model, the aver- 
age parameter likelihood over the region bounded by the 
fi^o and H to constraints, with a cut-off at A = 0.01, 
is £ avc ~ 0.92£j? for the full P98 data set. The aver- 
age parameter likelihood in the small A region, where 
the potential is effectively a cosmological constant model 
whose value depends on Va only, is C avc ~ 1.5£r. By 
a similar argument to above, the Bayes factor for the 
double-exponential potential model versus the PNGB 
model varies slightly over a range 0.92 < B < 1.5 for 
the full P98 data set. Similarly, 0.38 < B < 1.5 for the 
reduced "fit C" data set. Again, the Bayes factor is too 
small to confidently discriminate between the two mod- 
els. 

For the inverse power-law potential model, as shown 
in Fig. g, the fi^o and Hoto contours remain nearly par- 
allel towards the upper right-hand region where both 
Va and a are large. The average parameter likelihood 
over the region bounded by the Q^o and Hoto con- 
straints with arbitrary cut-offs at k" +2 Va/Hq = 10 15 
and a = 15 is £ avc — 0.14£i{ for the full P98 data set, 
and jCavo — 0.11£_rc for the reduced "fit C'data set, 
where Crc is the average parameter likelihood for the 
PNGB model based on the "fit C" data set. This gives 
a Bayes factor for the PNGB model versus the inverse 
power-law potential model of B = 7.2 for the full P98 
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data set, or B — 9.2 for the "fit C" data set. Both val- 
ues are large enough to provide slightly positive evidence 
that the P98 data set favours the PNGB model over the 
inverse power-law potential model. Since the Bayes fac- 
tor is only weakly positive, however, we cannot not place 
strong confidence in this conclusion. Note that the fl^o 
and -ffo^o contours both extend towards the smaller like- 
lihoods region. Therefore, increasing the prior ranges 
for Va and a will decrease the average parameter likeli- 
hood, and increase the slightly positive preference for the 
PNGB potential as compared to the inverse power law 
potential in the Bayes factor test. 



double-exponential potential model. 

For the inverse power-law potential model, the av- 
erage parameter likelihood over the region bounded by 
the f2^o and -ffo^o constraints with arbitrary cut-offs at 
K a+2 V A /Hg = 10 15 and a = 15 is £ avc ~ 1.9 x lCT 10 /^- 
This gives a Bayes factor for the PNGB model versus the 
inverse power-law potential model of B ~ 10 6 . This pro- 
vides very strong evidence that data set A would favour 
the PNGB model over the inverse power-law potential 
model. 



IV. DISCUSSION 



B. SNAP data set 

We have shown that the P98 data set does not particu- 
larly favour any of the PNGB, simple exponential poten- 
tial, and double-exponential potential models relative to 
the others. The P98 data set does disfavour the inverse 
power-law potential model as compared with the PNGB 
model. This might be considered to result from the fact 
that the 2cr-confidence region only overlaps with a small 
region (the low a region) of the whole parameter space 
allowed by the priors set on fi^o and ifo^o- 

In this section we want to compare the typical results 
we should expect when a future SNAP data set is used. 
We will use data set A, simulated from a fiducial PNGB 
model assuming no luminosity evolution. We will study 
the potential of this data set to distinguish the PNGB 
model from other quintessence models. 

For the simple exponential potential model, the aver- 
age parameter likelihood over the region bounded by the 
n,po and -ffo^o constraints, with a cut-off at A = 0.1, 
is £ avc ~ 2.1 x lCr 5 £ , where C = exp(-1910/2). 
In the small A region, we estimate the average param- 
eter likelihood to be £ avo — 6.6 x 10~ 31 £o- Compar- 
ing with the average parameter likelihood for the PNGB 
model, £ avc ~ 7.1 x 10~ 4 £o, the Bayes factor for the 
PNGB model versus the simple exponential potential 
model takes values in the range 34 - 10 27 , depending on 
the prior range for A. This would provide a very strong 
conclusion that data set A favours the PNGB model over 
the simple exponential potential model. 

For the double-exponential potential model, the aver- 
age parameter likelihood over the region bounded by the 
flc/,0 and Hoto constraints, with a cut-off at A — 0.01, is 
£ ave ~ 1.6 x 10 _3 £o. The average parameter likelihood 
in the small A region, is £ avo — 2.1 x 10~ 29 £o- The 
Bayes factor for the PNGB model versus the double- 
exponential potential model lies in the range 0.44 - 10 25 , 
depending on the the prior range for A. The lower bound 
of the Bayes factor, B = 0.44, is too small to confidently 
discriminate between the two models. However, provided 
small values of the parameter A are included in the prior 
range, the Bayes factor would provide a strong conclu- 
sion that data set A favours the PNGB model over the 



Let us now consider our results in relation to previous 
work concerning the feasibility of determining the prop- 
erties of the dark energy from future supernovae surveys 
p3,E3-pq|. One approach in past studies of scalar field 
quintessence has been to assume a potential V(<f>) over 
which the scalar field, <f>, would have slowly varied dur- 
ing cosmological time scales, and then test the efficacy of 
reconstructing the potential p2] , p4| . Another approach 
has been to assume that the quintessence field can be 
described by a perfect fluid with slowly varying equation 
of state P = w(z)p, expand w(z) in a power series, and 
then test the efficacy of determining the coefficients in 
the power series pa p8| . 

The conclusions of investigations to date are mixed 
p5| p8[ . Weller and Albrecht |2q] find that many mod- 
els can be distinguished with a fit to a linear equa- 
tion of state for the dark energy, P = w{z)p with 
w{z) = wq + w\z, but only if the current mass density, 
f2 m , is known to a high precision. Barger and Marfatia 
p7fl find that, even by putting Q m = 0.3 exactly, there 
is still a possibility of obtaining data sets which might 
not discriminate between quintessence and "fc-essence" , 
namely an alternative form of dark energy with a scalar 
field characterized by non-linear kinetic terms [ p7[ . Wang 
and Garnavich J2£| consider two classes of functions w(z), 
corresponding respectively to a linear variation and to k- 
essence. Using somewhat different techniques to other 
authors they are more optimistic about prospects for de- 
termining w(z) from future SNe la data. 

In the present paper we have taken a different approach 
from those above, by considering a class of models - the 
PNGB models - which are very well motivated from the 
point of view of particle physics, but for which the above 
methods will not always be satisfactory in the case of 
all plausible parameter ranges, given the potentially os- 
cillatory nature of w(z) and the corresponding fact that 
the scalar field may have varied over a wide range of val- 
ues of V{4>) over observable time scales. We fitted the 
simulated SNAP data sets to the exact di, (z) of different 
models obtained by numerical integration, and compared 
them to other models. 

One real drawback of all approaches is that as yet there 
is no preferred physical model for the dark energy. On 
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one hand this means that any approximations made in 
potential reconstruction methods may be too restrictive, 
since many different potential energy functions V(<j>) are 
conceivable, and many of these may give results degen- 
erate with each other. On the other hand, using a given 
Lagrangian for the quintessence field, as we have done, 
limits us to a model by model test. 

Nevertheless, we find that data sets such as those that 
would be produced by SNAP promise to be very suc- 
cessful on some tests, even if they will probably be less 
successful on others. In particular, while existing data is 
not yet sufficiently large to discriminate between various 
quintessence models or models with evolution [p|-|I(| , we 
have shown that the much larger size and smaller error 
bars of the simulated SNAP type IA supernovae data 
sets provide much tighter constraints on the parameters 
for quintessence models such as those corresponding to 
pseudo Nambu-Goldstone bosons. 

By evaluating Bayes factors in the context of the 
PNGB model, we have shown that future satellite SNe 
la data sets should have greater success in detecting 
whether the observed luminosity distance - redshift re- 
lation is purely cosmological in origin, or is significantly 
contaminated by evolutionary effects of the sources. The 
results of section IIB showed that although it may be dif- 
ficult to completely rule out luminosity evolution, if the 
true data were from a population with luminosity evolu- 
tion then this would provide a strong distinctive signal. 
We have only studied one simple illustrative supernova 
luminosity evolution function, but we expect that sim- 
ilar conclusions would apply to other simple luminosity 
evolution models. 

We have further shown that with the future data it 
should be possible to discriminate PNGB model from 
some other particular types of quintessence; in particu- 
lar, it gives a very different signature to that of simple 
inverse power-law potentials or simple exponential po- 
tentials. The case of a double-exponential potential gives 
a lower Bayes factor, and may therefore be more difficult 
to distinguish from the PNGB model. However, even 
in this case some distinction between the two models is 
possible if one allows a suitably large prior range of the 
parameter A, to include small values. 

A number of obvious extensions of our analysis are 
possible. In the case of testing source evolution versus 
the case with no evolution, for example, it would be in- 
teresting to determine by how much we can reduce the 
parameter j3 for the simulated SNAP data set B while 
still obtaining a very strong result for the discriminatory 
power of the relevant Bayes factor. Given the magnitude 
of the value obtained in section IIB, we suspect that the 
fiducial value of (3 could be significantly reduced. Like- 
wise many other tests could be performed with fiducial 
data sets based on other quintessence potentials. 

In conclusion, we find that future supernovae measure- 
ments such as those that would be afforded by the SNAP 
satellite, will have the power to significantly increase our 
knowledge of the properties of the dark energy in the 



universe. To be completely confident, however, we will 
require a deeper theoretical understanding of the nature 
of the dark energy and hopefully new input from funda- 
mental physics. 
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